Quantum Search by Local Adiabatic Evolution 
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The adiabatic theorem has been recently used to design quantum algorithms of a new kind, where 
the quantum computer evolves slowly enough so that it remains near its instantaneous ground state 
which tends to the solution pj|. We apply this time-dependent Hamiltonian approach to the Grover's 
problem, i. e., searching a marked item in an unstructured database. We find that, by adjusting the 
evolution rate of the Hamiltonian so as to keep the evolution adiabatic on each infinitesimal time 
interval, the total running time is of order vN , where N is the number of items in the database. 
We thus recover the advantage of Grover's standard algorithm as compared to a classical search, 
scaling as N. This is in contrast with the constant-rate adiabatic approach developed in |l|, where 
the requirement of adiabaticity is expressed only globally, resulting in a time of order N. 

PACS numbers: 03.67.Lx, 89.70.+C 



Introduction 

Although quantum computation is mostly a theoretical 
concept today, several quantum algorithms have already 
been designed and shown to outperform all known classi- 
cal algorithms, thereby giving a strong motivation to the 
development of quantum computers. Probably the most 
spectacular result is Shor's algorithm M, which can fac- 
tor a large number with a computation time polynomial 
in the size of the number, whereas all known classical 
algorithms require a (sub-)exponential time. Another re- 
markable algorithm, due to Grover, concerns the prob- 
lem of searching in an unsorted database |3| . Suppose we 
have a database of N items, one of which is marked. The 
goal is to find this unknown marked item by accessing the 
database a minimum number of times. Classically, N/2 
items must be tested, on average, before finding the right 
one. Grover's quantum algorithm performs the same task 
with a complexity of order y/N, giving rise to a quadratic 
speed up. 

While Grover's algorithm was originally presented 
within the standard paradigm for quantum computation, 
that is using a discrete sequence of unitary logic gates, we 
will turn here to another type of quantum computation 
where the state of the quantum register evolves contin- 
uously under the influence of some driving Hamiltonian. 
This concept of quantum computation viewed as a con- 
tinuous time evolution was pioneered by E. Farhi and S. 
Gutmann in [|4J, where they proposed an "analog" ver- 
sion of Grover's algorithm based on a time-independent 
Hamiltonian applied for a time T. Their algorithm re- 
quired T to be of order vN, which gives thus the same 
complexity as Grover's algorithm. In a more recent arti- 
cle with other coworkers, they considered an alternative 
class of continuous-time quantum algorithms based on 
a time-dependent Hamiltonian that effects an adiabatic 
evolution of the quantum register p|. If the evolution 
of a quantum system is governed by a Hamiltonian that 
varies slowly enough, this system will stay near its in- 



stantaneous ground state. This adiabatic evolution can 
be used to switch gradually from an initial Hamiltonian, 
whose ground state is known, to a final Hamiltonian, 
whose ground state encodes the unknown solution. The 
time required for this switching to remain globally adi- 
abatic determines the computation time. Farhi et al. 
have solved Grover's search problem using this adiabatic 
evolution approach, but this unfortunately resulted in a 
complexity of order N, that is no better than a classical 
algorithm that checks all possible solutions 0. 

In the present article, we show that one can recover the 
quadratic speed-up of Grover's original algorithm by con- 
tinuously adjusting the rate with which the initial Hamil- 
tonian is switched to the final Hamiltonian so as to fulfill 
the condition of adiabaticity locally, i. e. at each time. In- 
terestingly, this local adiabatic evolution approach makes 
it thus possible to improve the scaling law of the complex- 
ity of the quantum search algorithm simply by varying 
the evolution rate of the driving Hamiltonian. This offers 
the perspective of speeding up more sophisticated adia- 
batic evolution algorithms, such as those applied to NP- 
complete problems B . It might also be used to design an 
adiabatic evolution version of the nested quantum search 
technique proposed in Q for solving structured problems. 



Adiabatic theorem 

Consider a quantum system in a state \i/j(t)), which 
evolves according to the Schrodinger equation 



i Jt m)) = H(t)\m) 



(i) 



where H(t) is the Hamiltonian of the system (we let 
h = 1). If this Hamiltonian is time-independent and 
the system is initially in its ground state, then it will re- 
main in this state. The Adiabatic Theorem (see e. g. Q|) 
expresses that if the Hamiltonian varies slowly enough, 
roughly speaking, the state of the system will stay close 



to the instantaneous ground state of the Hamiltonian at 
each time t. More specifically, let \E k ;t) be the eigen- 
states of H(t), satisfying 



H(t)\E k -t) = E k {t)\E k -t) 



(2) 



where E k (t) are the corresponding eigenvalues and k la- 
bels the eigenstates (fc = labels the ground state). We 
define the minimum gap between the lowest two eigen- 
values as 



min \Ei (t) 

0<t<T 



Eo(t)) 



(3) 



and the matrix element of dH/dt between the two corre- 
sponding eigenstates as 



, dH . , . dH , , 

{—)x, = {E x] t\ — \E Q -t) 



(4) 



The Adiabatic Theorem states that if we prepare the sys- 
tem at time t = in its ground state \E ;0) and let it 
evolve under the Hamiltonian H(t), then 



|<£o;TbMT)}| 2 >l- £ 2 



provided that 



\/dH\ 
|\ dt i 



1,0 



<£ 



(5) 



(0) 



where e « 1. In particular, this implies that the mini- 
mum gap cannot be lower than a certain value if we re- 
quire the state at time t to differ from the instantaneous 
ground state by a negligible amount (a smaller gap im- 
plies a higher transition probability to the first excited 
state). This result can be used to design a new type of 
quantum algorithm based on a time-dependent Hamilto- 
nian |l|]. Assume we can build a Hamiltonian for which 
we know that the ground state encodes the solution of a 
problem. Then, it suffices to prepare the system in the 
ground state of another Hamiltonian, easy to build, and 
change progressively this Hamiltonian into the other one 
in order to get, after measurement, the sought solution 
with large probability. The Adiabatic Theorem imposes 
the minimum time it takes for this switching to be adia- 
batic, and this time can be thought of as the algorithm 
complexity. 



Global vs local adiabatic evolution for solving the 
quantum search problem 

We now apply this adiabatic evolution method to the 
problem of finding an item in an unsorted database. Be- 
fore describing our local adiabatic search algorithm, we 
first summarize the method of M , based on a global adia- 
batic evolution. Consider a set of N items among which 
one is marked, the goal being to find it in a minimum 
time. We use n qubits to label the items, so that the 
Hilbert space is of dimension N = 2™. In this space, 
the basis states are written |i), with i = 0, . . . , N — 1, 



while the marked state is denoted by \m). As we do not 
know \m) a priori, we use as an initial state an equal 
superposition of all basis states: 



N—l 



-y 



|Vt>> : 

We define two Hamiltonians: 

#o = J-hfoXVbl 



H„ 



\m)(m\ 



(7) 



(8) 
(9) 



whose ground states are \i/jq) and |m), respectively, each 
with eigenvalue 0. The time-dependent Hamiltonian un- 
derlying the algorithm is a linear interpolation between 
these Hamiltonians, that is 



H(t) = (l-t/T)H + t/TH„ 

or, with s = t/T 

H(s) = (l-s)H + sH m 



(10) 



(11) 



The algorithm consists in preparing the system in the 
state IV'(O)) = \ipo) arK l then applying the Hamiltonian 
H(t) during a time T. First, we notice that 



,dH _ ds dH 1 dH 

)i,o — -77 \— ;— )i,0 — 7^ (— ;— )i.o 



dt dt ds 1 ds 

Using Eq. (|l(j) or (^lj), we can solve the eigenproblem 
(0) and evaluate (H) and (H). The eigenvalues of H(s) 
are plotted as a function of s in Fig. [y. The highest 
eigenvalue E 2 = 1 is N — 2 times degenerated, while the 
two lowest ones Eq and E\ have a degeneracy one. The 
difference between these lowest two eigenvalues defines 
the gap g. We find 



AT-1 , 



N 



dH 

\-r- )i,o 
ds 



< 1 



(13) 
(14) 




FIG. 1. Eigenvalues of the time-dependent Hamiltonian 
H(s) as a function of the reduced time s, for N = 64. 



We see that the minimum gap g m in = 1/v A" is at- 
tained for s = 1/2, so that the adiabatic condition flq) is 
verified provided that 



N 
T> — 

e 



(15) 



Thus, the computation time is of order N, and there is no 
advantage of this method compared to a classical search. 
Now, let us see how to improve on this adiabatic evo- 
lution method. One should note that by applying Eq. 
(pi) globally, i. e. to the entire time interval T, we impose 
a limit on the evolution rate during the whole computa- 
tion while this limit is only severe for the times around 
s = 1/2, where the gap g is minimum. Thus, by divid- 
ing T into infinitesimal time intervals dt and applying 
the adiabaticity condition locally to each of these inter- 
vals, we can vary the evolution rate continuously in time, 
thereby speeding up the computation. In other words, we 
do not use a linear evolution function s(t) any more, but 
we adapt the evolution rate ds/dt to the local adiabatic- 
ity condition. Let us find the optimum s(t) with the 
boundary conditions s(0) = and s(T) — 1. Applying 
Eq. (ra) to each infinitesimal time interval, we get the 
new condition 



<£ 



9 2 (t) 



» ds > 



1,0 



(16) 



for all times t. Using Eqs. ( |l3| ) and (|14|), we choose to 
make the Hamiltonian evolve at a rate that is solution of 



d i = £g2(t) = £ 



N- 1 , 
1-4 s(l-s) 

N V ' 



(17) 



where e << 1. After integration, we find 

1 N 



t 



2e VAT- 1 



arctan 



(\A/V-1(2.< 



arctan \/N — 1 



(18) 



By inverting this function, we obtain s(t) as plotted in 
Fig. g, which shows the gradual change in the switching 
between Ho and H m . We see that H{t) is changing faster 
when the gap g(s) is large while it evolves slower when 
s is close to 1/2, that is where the gap is minimum. We 
may now evaluate the computation time of our new algo- 
rithm by taking s = 1 . With the approximation N » 1, 
we obtain 



T^VN 



(W) 



improving upon Eq. ([Tq). As a consequence, we thus 
have a quadratic speed-up compared to a classical search, 
and this algorithm can be viewed as the adiabatic evo- 
lution version of Grover's algorithm. In the Appendix, 
we show that this algorithm is optimal, that is the com- 
putation time cannot be shorter than 0(yN~) using any 
other evolution function s(t). 




FIG. 2. Dynamic evolution of the Hamiltonian that drives 
the initial ground state to the solution state: plot of the evolu- 
tion function s(t) for N = 64. The global adiabatic evolution 
method of |l| would appear here as a straight line between 
s(0) = and s(l) = 1. 



Conclusion 

In this article, we have applied the adiabatic evolution 
technique of H to design a quantum algorithm for solving 
Grover's problem, i. e., the search for a marked item in 
an unstructured database. We have shown that applying 
the Adiabatic Theorem globally (as in ]]]]) imposes a run- 
ning time of order N, where N is the number of items in 
the database, whereas adjusting the evolution rate of the 
Hamiltonian continuously in time so as to fulfill the adi- 
abaticity condition locally results in a time of order y~N. 
We therefore recover the advantage of Grover's usual al- 
gorithm compared to a classical search M. We should 
notice that this speed up was achieved by switching the 
Hamiltonian according to Eq. (|lj) , which is only possible 
because here the gap g(s) can be derived analytically and 
does not depend on the solution of the problem (in our 
case, \m)). As long as these conditions are satisfied, such 
a local adiabatic evolution method could be applied to 
more complicated - and more realistic - problems such 
as NP-complete problems, treated using either a quan- 
tum adiabatic evolution algorithm |3| or a nested version 
of Grover's algorithm exploiting the problem structure 
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Appendix: Proof of optimality 

In this Appendix, we show that using our algorithm, 
no other choice of the evolution function s(t) could lead 
to a better complexity than Eq. ( J19J ) . This proof follows 
closely the lines of the optimality proof of the "analog" 



Grover's algorithm based on a time-independent Hamil- 
tonian m. 

Let \ipmi t) be the state of our quantum register during 
the computation when the solution state is \m). After 
the computation time T, the states corresponding to dif- 
ferent solutions (m and mf) must be sufficiently different: 

l-|(V> m ,T|Vw,T)| 2 >e Vra^m' (20) 

Let us decompose the Hamiltonian H(s) in two parts: 

H(s) = #i(s) + H 2m (s) (21) 

where 

#r(s) = I-(l-s)\^ )(H (22) 

H 2m {s) = -s\m)(m\ (23) 

\ipm) aud \ipm') are solution of the Schrodinger equations 



then leads to 



i—\ip m ,t) = (Hx+H 2m )\ip m ,t) 

i—\lp m ',t) = {Hi +H 2 m')\^m',t) 

with initial conditions 



(24) 
(25) 

(26) 



We will now derive a necessary condition on T for Eq. 
( pcf ) to be satisfied. Using Eqs. ( j24| ) and (J2^), we have 



-[1-|(VW|Vw,*)| 2 ] 



(27) 



= 2Im[(Vw*l#2m-#2m#m^)(VW,#m,*)] ( 28 ) 

< 2\(ip m ,t\H 2m -H 2m ,\iP m ,,t)\\(^ m , 1 t\^ m ,t)\ (29) 

< 2 [|(VW,t|#2m|VW,*)| + |<VWl#2m'|VW,t>|] (30) 

Summing over m and m', we get: 



- ^2[i-\{ip m ,t\i> m >,t)\ 



dt 



m,rn' 



< 



4 Yl \{^m,t\H2m\lpm',t)\ 



m,m' 



< 4J2 ll#2 ro |^ m ,t>|||||</W,t}|| 



m,m' 



< 4/v£||ff 2m |Vm,*)|| 



(31) 
(32) 
(33) 
(34) 



where we used the Cauchy- Schwartz inequality along 
with the fact that \ip m i,t) is normalized. The property 

J2 \\H 2m \ib,t)\\ 2 = s 2 =► ]T \\H 2m \^t)\\ < VNs (35) 



- ^[1H(VW|VW)| 2 ]<4A^ S 



dt 



(36) 



We may now integrate this inequality using the initial 
conditions (E6f): 



J2 [I" IW< m ,*|Vw,i)| 2 ] < 4AViV / s(t)cft (37) 



Finally, using condition (20) and < s(t) < 1, we find 

(38) 



£JV- 1 

4 JN 



or, with A >> 1, 

T> -y/N 
~ 4 



(39) 



We conclude that, in order to be able to distinguish be- 
tween states corresponding to different solutions (e > 0), 
the computation must last a minimum time of order vJV, 
which is what we found in Eq. (|19|). Our choice of s(t) 
is thus optimal. 
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